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The nonlinear redshift space probability distribution 
function in models with local primordial non-Gaussianity 
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ABSTRACT 

We use the ellipsoidal collapse approximation to investigate the nonlinear redshift 
space evolution of the density field with primordial non-Gaussianity of the local /re- 
type. We utilize the joint distribution of eigenvalues of the initial non-Gaussian shear 
field and evaluate the evolved redshift space probability distribution function (PDF). 
It is shown that, similar to the real space analysis, the underdense tail of the nonlinear 
redshift space PDF differs significantly from that for Gaussian initial conditions. We 
also derive the lowest order correction of the Kaiser's formula in the presence of a 
non-zero /„;. 
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1 INTRODUCTION 

Cosmological probes of primordial non-Gaussianity have 
recently attracted much attention because of their potential 
ability to discriminate between d i fferent inflationary mod- 
els (e.g.. iBuchbinder etafl |200S| ; iKhourv fc Piazza! l20Qg| ; 
ISilvestri fc Trodden! 2008L and references therein). Con- 
straints on primor dial non-Gauss i anity mainly come 
from the CMB jKomatsu et all [20091; Ifflkage et all 



non-Gaussian model, in which the primordial perturbation 
potential is 



* = <l> + f«l(<P -(<?))■ 
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This paper is concerned with one particular measure of 
large scale structures: the probability that a cell of volume 
V, placed at random in the nonlinear redshift space density 
field, contains a certain amount of mass (or, equivalently, 
is denser than the background by a certain amount). This 
statistic is known as the nonlinear redshift space probability 
distribution function (PDF). Our goal is to estimate this 
distribution for scales as small as a few Mpc in the local 
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Here, $ is the Bardeen potential, 4> is a Gaussian poten- 
tial field and /„; is the nonlinear quadratic parameter. The 
right-hand side of eq. {TJ shows the first two terms of an (in- 
finite) Taylor series in (j>. However, since \<l>\ ~ 10 -5 , one usu- 
ally ignores higher order corrections, and commonly refers 
to this simplified model as the f n i model. This definition of 
$ is consistent with most of the rece nt studies on the lo- 
cal f n i model (b ut our earlier studies, lLam fc Sh cth 2009; 
lLam et al.ll2009l . defined "1? as the Newtonian potential). 

Our approach is based on previous work which devel- 
ops the formalism needed for estimating the evolution of the 
density PDF fr om Gaussian initial conditions in real and 
redshift space l|Lam fc Shethl l2008al lbl). The evolved PDF 
depends on the collapse dynamics and the statistical proper - 
ties of the initial density field. Recentlv. lLam fc Shethl (2009) 
used the fact that only the initial conditions are affected by 
primordial non-Gaussianity to model the evolution of the 
real space nonlinear PDF for the local non-Gaussian model. 
Their approach provided good quantitative agreement with 
measurements in numerical simlations. In what follows, we 
will assess whether this is also true in redshift space. 

Although it is possible to study the evolution 
of the redshift space PD F usin g perturbatio n the- 
ory methods jBernardeau Il994 ; Hivon et al.l Il995l ; 

IScoccimarro fc Friemanl Il999l ; \ Bernardeau et ahl |2002| ). 
this has, somewhat surprisingly not been extended to the 
local non-Gaussian model. Thus, it is not obvious how 
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Kaiser's formula relating the variance (the second mom ent 
of the PDF) in real and redshift space l|Kaiserl 119871 ) is 
modified when the initial conditions are non-Gaussian. 
(Although Kaiser's original derivation makes no explicit 
assumption about Gaussianity, the Gaussian assumption 
play s an important role in other de rivations of his formula, 
e.g., iFisherl 1 19951 : lOhta et alj|2004 l Our approach is quite 
different from Kaiser's, as it is based on an approximate 
model for the dynamics - the ellipsoidal collapse model 
- which reduces to perturbation theory at early times 
l|Bond fc Mversl Il996l). but allows one to study mor e 
nonlinear structures llSheth et al. I l200ll : |Pesiacquesll2008ft . 
lLam fc Shethl (|2008bh showed that the ability to probe 
deeper into the nonlinear regime, using a dynamical model 
that does not assume spherical symmetry, was crucial for 
modelling the PDF, especially in redshift space. However, 
implementing this approach requires knowledge of the 
initial shear field. For Gauss ian initial condition s, this has 
been known for some time l|Doroshkevichlll970l ). but how 
Doroshkevich's formulae are modified for the local f n i 
model has been shown only recently (|Lam et al.l [2009). 
Hence, we now have the necessary ingredients to study the 
redshift space PDF. 

Properties of the initial shear field in the local non- 
Gaussian model are briefly reviewed in Section 12.11 The 
dynamics of ellipsoidal collapse and the calculation of the 
nonlinear redshift space PDF are described in Sections 12.21 
and 12.31 We compare our model predictions with numeri- 
cal simulations in Section [3] We summarize our results in 
Section 2] A perturbative treatment of our model is given 
in Appendix [X] this shows explicitly that Kaiser's formula 
holds to lowest order in /„;, but that at higher order, it is 
modified. 



2 THE REDSHIFT SPACE DENSITY PDF IN 
THE LOCAL NON-GAUSSIAN MODEL 

Let us define the nonlinear overdensity of a region of volume 
V containing mass M by 

M_ 

W 

where p is the mean density. We will use p 3 to denote the 
corresponding quantity in redshift (rather than real) space. 
This section studies the expected dependence of the PDF 
of p s on the value of f n i, when the primordial potential is 
given by equation (JTJ). 

To proceed, we use the assumptions made when dealing 
with Gaussian initial distributions (/„; = 0): there is a local 
mapping from the eigenvalues \j of the initial deformation 
tensor to the nonlinear overdensity p s (see section |2"1?)) : and 
statistics on the smoothing scale V at the present time are 
related to statistics on a different smoothing scale in the 
initial conditions - the relevant initial smoothing scale is 
the one which contains the same mass (so it is larger for 
overdense cells, and smaller for underdense cells). 

Therefore, the nonlinear redshift space PDF of p s is 
given by 



1 + 8 



(2) 



p s p{p 3 \V)dp 3 



dAdep(Ajcr) <5 D [p s = p' s (\ e) 



(3) 



size V), A denotes the 3 eigenvalues (our convention is to 
have Ai > A2 > A3) of the initial 3x3 deformation tensor 
when smoothed on scale M (not V), a 2 denotes the vari- 
ance of the initial density fluctuation field on this smooth- 
ing scale (the initial density fluctuation Si is defined by 
Si = TrA), p' s (X,e) is the local mapping from the initial 
field to the evolved density given by the ellipsodial collapse 
model (spherical evolution models assume that the mapping 
is driven by the initial density Si only) , and e represents the 
rotation vector from the line-of-sight direction to the prin- 
ciple axis of the el lipsoid. Equation (|3|l h as the same form 
as equation (8) of lLam fc She th (2008bj) but, m our case, 
p(X\a) is the joint distribution of the initial eigenvalues Aj 
in the f„i model rather than in the Gaussian model. 

Before we compute p(A|<r), note that equation © does 
not guarantee a properly normalized PDF. To ensure the 
correct normalization, we set p' — Np and p' 2 p(p') = p 2 p{p) 
where N is chosen so that both f dp' p(p') and f dp' p'p(p') 
equal unity (|Lam fc Shethll2008ah . 



2.1 Initial conditions in the f„i model 

Let p(X\8i,a) denote the distribution of the Xj at fixed 81, 
and let po(A|cr) and po(X\Si, a) denote the corresponding 
quantities when f n i = 0, i.e., for Gaussian initial conditions. 
(Note th at this mean s pp(Si ) is a Gaussian.) One of the main 
results of lLam et al. I (|2009h was to show that 



p{X\8i,a) 



po{X\8i,(j) 
3 4 /4 



5 
2a 2 



5/2 



exp 



5Sf 151 
2a 2 + 2a 2 



T(5/2) 

x (A1-A2XA2-A3XA1-A3), (4) 

where the final expression for po is from iDoroshkevichl 
(1970), and I is the sum of the three permutations of XiXj 
where i j= j. Therefore the joint distribution of Xj in the /„; 
model is 



p(X\a) = p(8i\a)p(X\Si, 



p(Si\o-)p (X\Si,a), (5) 



where p(<5;|cr) is the distribution of the linear overdensity in 
the /„; model, and po(X\5i,a) is really a function of Xi/a 
and 5i/a. For the /„; values of current interest, p(<5;|<r) is 
only weakly non-Gaussian, so it can be approxim ated by 
the Edgeworth expansion (e.g.. lLam fc Shethl2009l . who also 
discuss the limitations of this approximation). Hence, the 
joint distribution of Xj is 



p(X\a) = 



1 + ^H 3 (Si/a] 



Po(A|ct) 



(6) 



where Hi{5i/a) = (81 /a) 3 — 3(5; /a) is the Hermite poly- 
nomial. The depend ence on f n i is encoded in the skewness 
parameter 0-5*3 (e.g.. IScoccimarro et al. | |2004 and note that 
our convention means that S3 is of same sign to f n i)- As a 
result, equation (J3j) becomes 



p s p(p s \V) dp s = / dAde 



1 + ^lH 3 (5i/a) 
b 



x po(Ajo-) 8 [p s = p' s (A, e)] 



where p 3 = M/M (M = pV is the average mass in cells of 



Except for the term in square brackets (the Edgeworth cor- 
rection factor) , the quantity in the integral is the same as in 
the Gaussian case. If we think of this extra factor as a weight, 
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then the resulting nonlinear redshift PDF in the f n i model 
is just a suitably weighted version of that in the Gaussian 
case. The weight depends on a S3, i.e. on /„;. 

We can gain some intuitive understanding of the ef- 
fect of a nonzero f n i as follows. For f n i = —100, a S3 < 
so that overdense regions are suppressed compared to the 
Gaussian case (the weight factor is less than unity), whereas 
underdense regions are enhanced compared to the Gaus- 
sian case (the weight factor is larger than unity). Finally, 
note that 0-53 <C 1; on Mpc scales a S3 ~ —0.03 for 
f n i — —100, and it is a weakly decreasing function of scale 



(e.g., [Scoccimarro et 
what follows. 



a weakly 
all 120041 ). 



This will be important in 



2.2 Ellipsoidal collapse and the nonlinear 
overdensity 

The next step is to estimate how p s depends on (A,e). In 
real space, the ellipsoidal evolution model sets 



3 n 

TT — 



(1-V3) : 



7=1 



Re (l-St/ScY 



(8) 



7=1 



l|Lam fc Shethl [2008a) , where Rj are the initial lengths of 
the patch which is now an sphere of radius Re, and S c is 
the critical value of spherical collapse model (its exact value 
depends weakly on cosmology: S c ~ 1.66 for the ACDM 
cosmology for which we show simulation data in the next 
section). In redshift space, the model sets 



81 

Pr 



E 



/ {Ri\ k - A* h Si [l - (1 - Si/Sc)*''*- 1 ] /3} 



^ Rl(l - X h ) - A' h [1 - S t /3 - (1 - St/SJ'o/*] 

(9) 

l|Lam fc Shethl l2008bh . where / = dlnD/dlna 
with D(t) the linear growth factor, (ei,e2,es) = 
(cos tjj sin 9, sin ip sin 9, cos 9), and A' h = 3/ ~}2j(Rk)~ , 
where the Rj are the initial axis lengths. 

2.3 Nonlinear PDF in redshift space 

Equations (|3j) and Q are the bases for the computation 
of the nonlinear redshift space PDF. The analysis simplifies 
considerably if we approximate 0S3 as a constant for a given 
Eulerian smoothing scale (recall that the scale dependence is 
rather weak): crS^ro, p s ) « crSziro). With this assumption 
we can write equation <[3j as 

p 3 p(p s \V)dp s = j d\dep(X\l) 5 D [p s = p' s (a\,e)] 

= JdXde l+^H 3 (8i) 

x p (A|l)fo [p s = p' s {a\,e)] . (10) 

In practice, we construct the PDF by Monte Carlo solution 
of the integral. This is straightforward because the six inde- 
pendent components of the deformation tensor <t>ij = didj<& 
can be combined in the form C = {x, y, z, $12, $23, $31 }, 
where 



V 



$11 - $2 



$11 + $2 



2$33 



(11) 



The reason for doing this is t hat, to second o rder in /„;, 
only x has non-zero skewness (|Lam et alj|2009h . Therefore, 
we can draw the other five parameters from Gaussian dis- 
tributions with variance 

(V 2 ) = = ^, and <* 2 } = y. (12) 

The parameter x is drawn from an Edgeworth distribution. 
The associated (Ai, A2, A3) can be computed by solving the 
eigenvalue problem and the nonlinear redshift PDF is then 
evaluated using equation (jTDJ . 

In the next section, we compare this full solution with 
measurements in simulations. Note however that, in the limit 
p s — 1 -C 1, the nonlinear redshift space PDF can be solved 
perturbatively. Appendix[X]provides details and shows that, 
to lowest order, the variance in the redshift space counts is 
related to that in real space by Kaiser's formula for /„; = 0; 
the dependence on /„,; enters at higher order. 



3 COMPARISON WITH SIMULATIONS 

We now compare the predictions of our model with measure- 
ments of the non l inear PDF in numerical simulations from 
iDesiacques et all |2009l ) . The numerical simulations followed 
the evolution of 1024 3 particles in a periodic cube of sides 
1600/i _1 Mpc. The background cosmology was ACDM with 
(n m ,Q b ,n a ,h,a 8 ) = (0.279,0.0462,0.96,0.7,0.81). 

Figure \T\ compares our model with the measured red- 
shift space PDF of counts in 8/i _1 Mpc spheres. In the up- 
l per panel, the solid symbols show the PDF for f n \ = 0. We 
have not shown results for f n \ = ±100 in that panel since 
they only slightly differ from the Gaussian case. Instead, 
the symbols in the bottom panel show the fractional devi- 
ation in these models relative to the Gaussian case. Filled 
and open symbols are for f„i = 100 and /„; = —100, re- 
spectively. As we can see, a positive f n i slightly skews the 
PDF towards overdense regions. Conversely, the fraction of 
underdense regions is enhanced for negative /„;. In this re- 
spect, the redshift space PDF shows the same qualitative 
dependence on /„; as the real space PDF, as expected (c.f. 
discussion following equation [TJl . 

The dashed, solid, and dotted curves in the top panel 
show the predictions for /„j = 100, 0, and —100 respectively. 
The differences are small, so the curves appear almost identi- 
cal, but the bottom panel shows that they are indeed slightly 
different from one another, and that our model provides a 
good description of the ratios, except in the high density tail 
where it underpredicts the dependence on /„; . Note however 
that, in this strongly nonlinear regime, our model drastically 
overpredicts the Gaussian counts. 

We believe we understand why our model is more suc- 
cessful at predicting the ratio than the counts themselves. 
This is because at least some of the discrepancy at p s 3> 1 
arises from the fact that the highly nonlinear virial mo- 
tions within halos will act to erase large density contrasts - 
these motions are not part of our model. Figures 1 and 2 in 
lLam fc Shethl l|2008ri ) show that, for f nt = 0, virial mo- 
tions reduce the p s ^> 1 tail, enhance the intermediate 
p s ~ 2 region of the PDF, and have almost no effect on 
the p s < regime. Removing virialized motions within ha- 
los from the measurements substantially reduces the dis- 
crepancy between theory and the simulations at p s > 0. 
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Figure 1. Comparisons of the measured PDF with our model for counts in 8h~ x Mpc sphere. The upper panel shows log(pp(p)) against 
ln(p) for the measured PDF (solid symbols, f n i = 0) and the theoretical prediction obtained by evaluating the nonlinear PDF eq. for 
fni = (black, solid), —100 (magenta, dotted), and 100 (magenta, dashed) respectively. The lower panel shows the logarithm of the ratio 
between the f n i and Gaussian counts. The filled and empty symbols indicate the measurement for f n i = 100 and —100, respectively. 
The predictions obtained by applying the Edgeworth expansion weighting are represented by the dot-dashed (red, f n i = —100) and the 
short-long-dashed (red, /„; = 100), respectively. 



Although virial motions do not depend on f„i, their net ef- 
fect depends upon the halo mass function. Since the later 
does depend on /„;, we may thus expect a slightly stronger 
suppression when f n i > (as the abundance of massive ha- 
los is slightly enhanced). On the other hand, the real space 
PDF, which is also affected by virial motions, has a more 
pronounced high density tail for /„; > 0. As a result, the ra- 
tio of non-Gaussian counts to Gaussian counts depends only 
weakly on /„;. Therefore, our model can provide a reason- 
able description of the ratio even though it fails at describing 
the high p s tail of the Gaussian density PDF. 



Figure [2] shows a similar comparison on smaller scales 
(spheres of radius 4/i _1 Mpc). The dependence on /„; is 
smaller compared to the previous figure. Our model still pro- 
vides a good description of the ratio relative to the /„; = 
counts, except at the highest densities where it overpredicts 
the f n i = counts and underpredicts the ratio. Note again, 
that we expect much of this discrepancy to be reduced if we 
were to remove virial motions from the simulations. 



4 DISCUSSION 

We used the ellipsoidal evolution model to study the redshift 
space probability distribution function of the nonlinear dark 
matter density field in the local non-Gaussian /„; model. 

A perturbative analysis of the density PDF eq.Q shows 
that, at the lowest order, Kaiser's formula still holds in the 
fni model (although his original derivation does not assume 
Gaussianity explicitly, other derivations of the formula have 
done so, as discussed in the Introduction). The effects of 
fni 7^ appear in the first order corrections to the variance 
(and higher order moments). One could, therefore, constrain 
fni from large scale structure by measuring the variance and 
the higher order moments and comparing with the pertur- 
bative quantities in Section [X] (with some dynamical models 
to determine Ui). 

Our approach remains accurate on smaller scales 
where perturbative treatments are not useful. Simulations 
show that the dependence on f n i is qualitatively similar 
to that for the real space PDF: for positive /„( (posi- 
tive a S3) both PDFs skew slightly towards overdense re- 
gions. In addition both show stronger /„; dependence in 
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Figure 2. Same as figure[T] but for counts in spheres of radius Ah ^^Mpc. 



the underdense regions, suggesting that void abundances 
should be good probes for primordial non-G aussianity (e.g., 
iKamionkowski et al.l 120091 ; lLam et all 120091 ). Our model 
(equation [7]) captures these trends (Figures [T] and [2} ■ Since 
it is explicitly a redshift space calculation, it would be inter- 
esting to see if it correctly predicts the /„; dependence of the 
PDF of the flux in the Ly-q forest, that has recently been 
simulated bvlViel et alJ (|2009h . This work also provides the 
foundation for constraining f„i in future galaxy surveys (e.g. 
the change in the redshift space halo/galaxy power spec- 
trum by_TOmbj^iinTwitli_th«_s£ale d 

bias ijDalal et al.ll2008l : ISlosadl2009l : lDesiacaues et alll2009h ). 
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APPENDIX A: PERTURB ATIVE TREATMENT 
OF REDSHIFT SPACE DISTORTIONS 

Equation © allows for a novel estimate of how redshift 
space statistics are expected to differ from those in real 
space as /„; varies. This is because the overdensity in red- 
shift space is 
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jLam fc Shetbll2008bl) . Here, fi = d\nD 1 /da w fi ' 55 where 
Di is the linear growth factor, fa — d\nD2/da where 
D 2 /Dl w -(3/7)n _1/143 , and v 2 » 1.62 and v z « 3.93 
are related to the spherical evolution model. 

Note that setting /„; 7^ simply changes the values of 
the averages over the As. Hence, to lowest order, 
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this is Kaiser's formula, so the relation between real and 
redshift space variance is unchanged from the Gaussian case. 

Of course, f n i matters for the higher order moments. 
The next higher order of the redshift-space variance is 
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The origin of these terms can be understood as fol- 
lows. When f n i = 0, then one can think of the three terms 
in Kaisers expression as being due to the density-density, 
density-velocity and velocity-velocity power spectra. Now, 
velocities are related to first derivatives of the potential, 
whereas densities are related to second derivatives. So one 
expects the lowest order corrections to the Gaussian result 
to scale as /„;. Terms in the first order correction (second 
equality in equation IA6[) can be interpreted as Bddd, Bdw, 
Bddv, Bddv, B vvv , and Bdw respectively, where B denotes 
bispectra, d and v axe density and velocity. Notice that this 
first order correction in the redshift variance is of lower order 
compared to the case where f n i — (and hence uSs = 0). 
This is generic for models w ith non-vanishing initial skew- 
ness ()Bernardeau et al1l2002l ) . 

This approach can be extended to estimate the real- 
redshift large scale relation in higher order statistics, 
for example the b ispectrum of galax i es to constrain f n i 
(see f or example, Scoccimarro et al.l (|2004l ) or more re- 
cently IJeong fc Komatsul ( 20091 ')'). However complications 
arise when one i nclud es scale depe n dent halo/galaxy 
bias dDalal et alj l20Qg|; ISlosaj 120091 ; iDesiacques et all 
120091 ; IDesiacques fc Shethl I2009T ) and the validity of the 
peak-background sp lit approach in computing halo bias 
llManera et al. 2009). These are beyond the scope of this 
paper and will be explored in future studies. 



